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Abstract 

We extend existing renormalization group calculations for the exponents describing 
scaling of star polymers and polymer networks constituted by chains of different 
species (the so-called copolymer star exponents). Our four loop results find appli- 
cation in the description of various phenomena involving self-avoiding and random 
walks that interact. 
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1 Introduction 



Star polymers, i.e. structures of polymer chains that are chemically linked 
with one end to a common core attain recently much attention both in theory, 
experiment, and technical applications [1]. On the experimental side this is 
due to their effective interaction properties that interpolate between those of 
colloidal particles and single chain polymers. For theory and simulation they 
present the most simple albeit non-trivial examples of polymeric networks 
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which, as shown by theory, provide all information to derive the scaling prop- 
erties of general polymer networks. For this reason the scaling properties of 
stars of long flexible polymer chains in a good solvent are of central interest 
[2,3]. The scahng of single long flexible polymer chains is perfectly described 
by a model of self-avoiding walks resulting in universal scaling exponents that 
only depend on the space dimension d [4] . However, the corresponding scaling 
exponents for star polymers will also depend on the number of arms / con- 
stituting the star [2]. Building a polymer star from chains of different species 
one finds an even richer scaling behavior [5,6]. The partition function of such 
a so-called copolymer star consisting of /i chains of species 1 and /2 chains of 
species 2 (each of the same size R) will show power law scaling in R according 
to [6]: 



with a family of copolymer star exponents 77/1/3 f*^^ /ii /2 = 0, 1, 2, . . .. In (1) £ 
is a microscopic scale and a fugacity factor is omitted. We consider a species 
A of self- and mutually avoiding chains (SAW), species B and B' of both self 
and mutually transparent chains (RW) and a species C of self-transparent 
but mutually avoiding chains (MAW). Copolymer stars are built from one 
or more of these species with mutual avoidance between any two species. 
Of special interest are the exponents of (a) homogeneous polymer stars rj^^, 
(b) symmetric and unsymmetric copolymer stars r/^^j , and ifj^^gi and (c) 
of stars of mutually avoiding walks Vfc^^- These exponents naturally arise 
when describing e.g. the scahng properties of the effective interaction of star 
polymers, the correlations of the vertices of a branched polymeric structure, 
the diffusion of particles near an absorbing polymer, and the transition of 
double stranded to single stranded DNA [1,7,8]. Moreover, the star exponents 
VfAfB uniquely define the scaling of copolymer networks of arbitrary but fixed 
topology [2,3,6]. The above mentioned facts together with theoretical reasons 
(e.g. the copolymer star exponents are related to the composite operators of 
polymer field theory and possess multifractal properties [9] which seemingly 
is in contradiction with stability conditions [10]) explain the high interest in 
their analysis. 

Currently, the field theoretical renormalization group (RG) [11] is recognized 
as the most reliable analytic tool to calculate numerical values of scaling ex- 
ponents in a systematic perturbative way. For the universality class of 0{m)- 
symmetric 0^ theory exponents are known up to high orders of perturbation 
theory: to the 5th order in s = 4 — d expansion [12] and to the 6th/7th order 
when calculated directly at fixed space dimension d = 3 [13]. However, the star 
exponents are known so far only to the 3rd order of perturbation theory [3,6]. 
A principal reason for this gap is that copolymer star exponents were intro- 
duced only recently [6] and that these exponents have no direct counterparts 
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in the theory of magnets. 

Our main message is given by the analytic expressions for the copolymer star 
exponents that we obtained by state-of-the-art RG methods to £^ accuracy in 
an ^-expansion. This is merely one order below the best known £-expansion for 
the critical exponents of the m- vector model [12]. When evaluated numerically 
[14,15] by appropriate resummation techniques [11] the expansions for the ex- 
ponents lead to results in nice correspondence with the available exact (2D) 
and MC data. This paper is organized as follows: the next, 2nd section intro- 
duces the main expressions relating the copolymer star exponents to polymer 
field theory, section 3 presents the e-expansions for the three cases discussed 
above: ^/^/^, , V^^fg ^fc^^' ^^^^ section we conclude and discuss 

our results presenting a numerical estimate drawn from our expansions. De- 
tails of our calculations along with numerical estimates of different quantities 
involving star exponents are the subject of a separate publication [14] . 



2 The renormalization group relations 

Here, we sketch the framework of polymer field theory suited to describe 
copolymer stars. The starting point is the Edwards model [4] generalized to 
describe a set of /i + /2 polymer chains in solution with self- and mutual- 
avoiding interactions Uab- Within this model, the configuration of each chain 
a = 1, . . . , fi + f2 is given by a path r"(s) in d-dimensional space, < s < Sa, 
Sa is the so-called Gaussian surface of the path r"(s). The Hamiltonian Ti/^/a 
of the model reads [3,6]: 




with densities Pa(r) = /o " ds5(r — r"(s)). In Eq. (2), the first term describes 
the chain connectivity whereas the second one stands for the excluded volume 
interaction. The partition function of a star of /i + /2 polymer chains 2^(*-^i-^2) 
is calculated as a functional integral (/ '^'[ra]) of the Boltzmann distribution 
with the Hamiltonian (2) over all possible configurations of the paths r"(s). 
A product of 5-functions serves as a constraint ensuring all polymer chains to 
have a common endpoint: 

Z(*f^f^\S^}^J P[rjexp{-^7i/,^,(r„,{5j)}'n'5(r„(0)-ri(0)).(3) 
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The mapping of the Edwards model to a 0^ field theory with one coupling of 
0(m)-symmetry in the limit m = is well established [4]. The polymer system 
we are currently interested in is the so-called ternary polymer solution, where 
polymer chains of two species (1 and 2) are present. It is described by an 
Edwards model (2) with three couplings [16]: 

Uab = {Ull, U22, Ui2 = U21} ■ (4) 

The partition function of the copolymer star (3) corresponds in the field theory 
to a vertex function with the insertion of a composite operator Hai^a of a special 
symmetry [6]. Ultraviolet divergences occur when the bare vertex functions 
are evaluated naively [11]. We apply the field theoretical RG approach [6] 
to remove the divergences and to study the scaling properties of the vertex 
functions of copolymer stars. Since the theory is renormalizable one can collect 
all divergences in the so-called renormalization factors Z and define a finite 
theory with renormalized parameters preserving the structure of the original 
one. The renormalized couplings gab and the renormalization factors Z depend 
on the parameter k which gives the scale of the external momenta in the 
renormalization procedure. This dependence is expressed by RG functions 
defined by the following relations: 



In Eq. (7), Z^f^f^ stands for the Z- factor renormalizing the vertex function 
of a composite operator, which corresponds to the copolymer star partition 
function (3). In the fixed points of the RG transformation, given by a common 
zero of the /3-functions (5), (6) the anomalous dimensions r^/^/j of the composite 
operators define the set of exponents for copolymer stars. Explicit expressions 
for the star exponents will be given in the next section. 



3 £^-expansions for the exponents 

The zeros of the /3-functions (5), (6) correspond to the fixed points of the 
theory. To clarify their physical significance we here refer to the three classes 
of species A, B, and C introduced above and denote nontrivial values by g^, 
gQ, and glj. The fixed points of (5) are Pii{g^) = for species A, Paa{0) = 
for species B and C. The non-trivial fixed points of (6) are /3i2(0, 0, ^fg.) ~ 




(5) 



(6) 



(7) 
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for the combination of species B, B' and for the interchain interaction of 
species C and finally f^uiOsi 0) 9u) ~ ^'^^ ^^e combination of species A and B 
[16]. Analyzing the fixed point structure, one is led basically to the following 
physically different situations: in the first case we combine species B and B' i.e. 
both sets of /i polymer chains of species B and of /2 polymer chains of species 
B' are non-selfinteracting RWs, but mutually avoid each other (we refer to 
this as to the case "G"). In the second case, "f/", the species A selfinteract (/i 
SAWs) and the species B do not (/2 RWs). In the third case, ''MAW", only 
species C of fc mutually avoiding RWs is present. In our calculations we use 
the dimensional regularization and minimal subtraction scheme [11]. Most of 
the star graphs could be derived from the four-point diagrams of the <^'^ theory. 
At the four loop level, however, we encountered a number of independent star 
graphs, which we evaluated following the R*-operation scheme [11]. Details 
of our calculations will be presented elsewhere [14]. We do not display the 
expressions for the /9-functions here although they are of special interest e.g. 
to study the crossover behaviour of the ternary polymer solution [16]. Our 
results for the copolymer star exponents (1) read: 



4h = + f'Mi + /2 - 3)^ - + /2 - 3) [3C(3) + /i 

+/2 - 3] ^ - /i /2 [70 C(5) - 54 C(4) - 56 C(3) - 132 A 

J 16 L 

-132 /2 + 60 C(3) /i + 60 C(3) /2 + 18 C(4) /i + 18 C(4) 
+30C(5) /i + 30C(5) /2 + 44/2 _ g/s + 44/2 _ 5/3 + qq 
-8 C(3) - 8 C(3) fi - 20 C(5) - 20 C(5) /| - 42 C(3) h h 

-30 C(5) h /2 - 16 /i fl - 16 fl h + 125] ^, 



(8) 



^/i/. - /i - 3/2)g + fmi - 91/2 - 33/1 + i9>fi + 42/1/2 

+42 C(3) h f2 + 25)^ - fl [712C(3) + 969/i + 2463/2 

-936C(3)/i - 2652C(3)/2 - A56f^ + 6Af^ - 1050 + 108/| 
-2290/1/2 + 224C(3)/i2 + 5A0C{S)f^ + 1188C(3)/i/2 + 50AfJ^ 

+492/2/2 - 577] -^ + fi [99200 C(5) - 34176 C(4) - 77136 C(3) 
-1 4096 L 

-136137 fl - 321075 /a + 128528 C(3) /i + 286200 C(3) /s 
+44928 C(4) fl + 127296 C(4) /a - 99840 C(5) fi - 311280 C(5) ^ 
+94376 f^ - 26752 /f + 216098 f^ + 2688 ff - 46872 /| 
+3240 f^ + 456730 A A - 59712 C(3) /^ + 8320 C(3) f^ 
-122904 C(3) /2 - 10752 C(4) /^ + 10368 C(3) /| - 25920 C(4) /^ 
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-8960 C(5) + 9600 C(5) f! - 36720 C(5) /| + 21600 C(5) fl 
-282504 C(3) A f2 - 57024 C(4) A A - 47280 C(5) /i A 
-212496 /i /I - 199768 /f ^ + 24192 A /| + 27624 ff h 
C(3) h fl + 59328 C(3) fl A + 79920 ((5) /i /| 

+80160 C(5) A^ A + 47808 fl fl + 65825] (9) 



^MAw ^ _!)£ + /(/_ i)(2/ - 5)- - /(/ - 1) 8C(3)/ - 20/ 



e 

16 



-19C(3) + 4/2 + 25j - + /(/- 1) [782 / + 222 C(3) 
+114 C(4) + 120 C(5) - 230 / C(3) - 48 / C(4) - 310 / C(5) 



-314 /2 + 42 /^ + 58 /2 C(3) + HO /^ C(5) - 647 



256 



(10) 



In (8)-(10), C,{x) is the Riemann zeta-function. The above expressions recover 
the results of Ref. [6] and describe the scahng of star-hke combinations of 
different species: two sets of RWs B, B' (8); /i SAWs and A RWs, A and B 
(9); / MAWs, C (10). The case of the homogeneous star exponent is included 
via r]^^ = r^y^o- 



4 Conclusions and outlook 



The expansions (8)- (10) are the starting point for the construction of analytic 
expressions for a number of physical quantities describing various phenomena 
(we mention some of them in the section 1): the exponents govern the short 
distance interaction of star polymers in colloidal solutions [1,15] as well as 
the reaction rate of diffusion controlled reactions with traps or reaction sites 
attached to polymer chains or to a star polymer [1,5]. One more recent appli- 
cation concern scaling of distribution of denaturated loops and unzipped end 
segments for a model of DNA denaturation [8] . Moreover, due to the special 
convexity properties of the spectrum of copolymer star exponents [9,10] they 
can be re-written in the multifractal formalism, in terms of a spectral function 
and Holder exponents [5,7,9]. Such a program is well outside the scope of this 
Letter and is the subject of separate publications [14,15]. We also hope that 
our study will motivate the analysis of these and similar problems by other 
tools. 

As a final note, we want to emphasize another particular feature of the spec- 
tra of exponents (8)-(10) distinguishing them from the usual ^-expansions ob- 
tained within 0^ theory. The field theoretic perturbative expansions are known 
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[11] to be asymptotic at best and one should apply a resummation procedure 
to calculate the numerical estimates from these series. Such a resummation 
allows, in particular, to restore the convergence of the series in e. However, the 
series (8)-(10) posses also a strong /-dependence. The standard resummation 
technique, leading to high-precision estimates in the field-theoretical scheme 
will give similar accuracy for low values of / only. To show this and to give 
an idea of the numerical estimations, based on e.g. Borel resummation refined 
by a conformal mapping [11], we display one particular result in Fig. 1. Here, 
we plot our resummed results for the star exponent — 1 + ^{rfjQ — frjw) 
[6] with the Flory exponent u — 0.588 in d — 3. This exponent governs the 
scaling of the homogeneous star partition function Z^*^'' with respect to the 
number of monomers in the chain: Z^*f^ ~ N^f-'^. We compare the re- 
summed £^ data with the results of MC simulations [17,18,19,20] and with 
the resummed three-loop pseudo-£ expansion obtained within the massive RG 
scheme at fixed d — 3 [21]. One observes a nice correlation within successive 
perturbative approximations and MC simulations for small / and an expected 
growing discrepancy when increasing the number of chains. A recent MC sim- 
ulation [20] has determined values for the exponents with high precision. One 
sees from the figure 1 that the usual resummation does not meet these values 
for high /. We beheve that a more sophisticated technique should also include 
the known large-/ behaviour. 
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Fig. 1. Homogeneous star exponent 7/ at d = 3. Open squares: forth order (e ), open 
triangles up: third order (s^), open triangles down [21]: third order RG calculated 
at fixed dimension (d = 3). Results of MC simulations are shown by filled symbols, 
triangles up: [17], triangles down: [18], diamonds: [19], discs: [20]. 
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